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We searched for a sidereal modulation in the MINOS far detector neutrino rate. Such a signal 
would be a consequence of Lorentz and CPT violation as described by the Standard- Model Extension 
framework. It also would be the first detection of a perturbative effect to conventional neutrino mass 
oscillations. We found no evidence for this sidereal signature and the upper limits placed on the 
magnitudes of the Lorentz and CPT violating coefficients describing the theory are an improvement 
by factors of 20 — 510 over the current best limits found using the MINOS near detector. 

PACS numbers: 11.30.Cp,14.60.Pq 



Neutrinos have provided many crucial insights into 
particle physics, including the existence of physics be- 
yond the minimal Standard Model (SM) with the detec- 
tion of neutrino oscillations [l], Because oscillations 
are interferometric in nature, they are sensitive to other 
indicators of new physics. Such indicators include po- 
tential small amplitude signals persisting to the current 
epoch whose origin is a fundamental theory that uni- 
fies quantum physics and gravity at the Planck scale, 
nip = 10^^ GeV. One promising category of Planck-scale 
signals is the violation of the Lorentz and CPT symme- 
tries that are central to the SM and General Relativity. 
The Standard Model Extension (SME) is the comprehen- 
sive effective field theory that describes Lorentz (LV) and 
CPT violation (CPTV) at attainable energies [|]. 

The SME predicts behaviors for neutrino flavor change 
that are different from conventional neutrino oscillation 
theory. The probability of flavor change in the SME de- 
pends on combinations of L, the distance traveled by the 
neutrino, and the product of distance and the neutrino 
energy, L x E^. For conventional oscillation theory the 
transition probability depends only on L/E^. The SME 
also predicts that the neutrino fiavor change probabil- 
ity depends on the angle between the direction of the 
neutrino and the LV/CPTV field in the sun-centered in- 
ertial frame in which the SME is formulated Ex- 
periments like MINOS whose neutrino beam is fixed 
on the Earth, are well-suited to search for this behavior, 
which would appear as a periodic variation in the de- 
tected neutrino rate as the beam swings around the field 
with the sidereal frequency w© = 27r/(23''56'"04.09053'*). 

MINOS has a near detector (ND) located 1 km from 
the neutrino beam source and a far detector (FD) lo- 
cated 735 km from the neutrino source. Because of their 
different baselines, the ND and FD are sensitive to two 
separate limits of the general SME formulated for the 
neutrino sector. The predicted SME effects for baselines 
less than 1 km are independent of neutrino mass , and 
both MINOS [6.] and LSND [J] reported searches for these 
effects. Recent theoretical work has shown that SME ef- 
fects are a perturbation to the dominant mass oscillations 
for neutrinos having the appropriate L/E^, to experience 
oscillations [8]. Since the probability for transitions due 
to LV increases with baseline, experiments with baselines 
greater than « 100 km are especially sensitive to LV and 
CPTV. The following analysis using MINOS FD data is 
the first search for perturbative LV and CPTV effects in 



admixture with neutrino oscillations. 

According to the SME, the transition probabil- 
ity for i>f^ — >■ i^T- transitions over long baselines is 

P^r — P^iT + P^vT , where pjjP is the conventional mass 
oscillation probability for transitions between two flavors 
and P'l^T is the perturbation due to LV and CPTV, with 
pI^JIpI!P « 1. In the SME, P^^^ is given by Q 

+ {P!^^Jrt.SmU^Tg, + {P^2])^^C0SLJg,Ts, (1) 

+ (P^'^ Sin 2wer® + (P^'j )rM cos 2u;qTq | , 

where L = 735 km is the distance from neutrino produc- 
tion in the NuMI beam to the MINOS FD [2, T© is the 
local sidereal time (LST) at neutrino detection, and the 
coefiicients iP^\^, {P^^Jr^, (PA^rt., (Pe^Un and 
{Pl^'')Tfj. contain the LV and CPTV information. These 
coefiicients depend on the SME coefficients that explic- 
itly describe LV and CPTV, (ol)^^ and (ci)|^, as well as 
the neutrino mass-squared splitting, Am^2 H- For two- 
fiavor transitions, only the real components of the (ol)"^ 
and {cl)"^ contribute to the transition probability. 

The magnitudes of the functions in eq. ([T|) depend on 
the direction of the neutrino propagation in a fixed coor- 
dinate system on the rotating Earth. The direction vec- 
tors are defined by the colatitude of the NuMI beam line 
X = (90°- latitude) = 42.17973347°, the beam zenith 
angle 6 — 86.7255° defined from the z-axis which points 
up toward the local zenith, and the beam azimuthal angle 
(j) — 203.909° measured counterclockwise from the a;-axis 
chosen to lie along the detector's long axis. 

This analysis selected data using standard MINOS 
beam quality requirements and data quality selections 0, 
0]. The neutrino events used must interact in the 4.0 kilo- 
ton FD fiducial volume fs?] and be charged-current (CC) 
in nature. The selection method described in Q allowed 
the identification of the outgoing muon in a CC interac- 
tion. As in Q, we focused on these events to maximize 
the disappearance signal. 

The data used come from the run periods listed in 
Table H] also shown are the number of protons incident 
on the neutrino production target (POT) for each period 
and the total number of events observed. To avoid biases, 
we performed the analysis blindly with the procedures de- 
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termined using only the Runs I and II data. The Run III 
data, comprising more than 50% of the total, were in- 
cluded only after finalizing the analysis procedures. 

TABLE I: Run Parameters 





Run Dates 


POT 


CC Events 


Run I 
Run II 
Run III 


May05 - Feb06 
Sep06 - Jul07 
Sep07 - Jun09 


1.24 X 10^" 
1.94 X 10^° 
3.88 X 10^° 


281 
453 
954 



We tagged each neutrino event with the time deter- 
mined by the Global Positioning System (GPS) receiver 
located at the FD site that reads out absolute Universal 
Coordinated Time (UTC) and is accurate to 200 ns flo| . 
The GPS time of the accelerator extraction magnet signal 
defined the time of each 10 /iS beam spill. We converted 
the time of each neutrino event and spill to local sidereal 
time T0 (LST) in standard ways [11[ . The uncertainty in 
the GPS time stamps introduced no significant system- 
atic error into the analysis |6j. We placed each detected 
CC event into a histogram that ranged from 0-1 in local 
sidereal phase (LSP), the LST of the event divided by the 
length of a sidereal day. We used the LSP for each spill 
to place the number of POT for that spill, whether or 
not there was a neutrino event associated with it, into a 
second histogram. By dividing these two histograms, we 
obtained the normalized neutrino event rate as a function 
of LSP, in which we searched for sidereal variations. 

Since the search for sidereal variations used a Fast 
Fourier Transform (FFT) al gori thm that works most ef- 
ficiently for A/" = 2" bins [l2| and eq. ([U only puts 
power into the Fourier terms wqTq and 2a;©r©, we chose 
A/" = 2^ = 16 bins to retain these harmonic terms while 
still providing sufficient resolution in sidereal phase to de- 
tect a sidereal signal. With this choice, each bin spanned 
0.063 in LSP or 1.5 hours in sidereal time. 

The statistical similarity of the event rates for all runs 
was tested by comparing the rate for run i in LSP phase 
bin j, Rij, with the weighted mean rate for that bin, Rj. 
The distribution of r = (Rij — Rj)/aij, where cr.y is the 
statistical uncertainty in Rij for all is Gaussian with 
f — and cr = 1, as expected for statistically consistent 
runs. Given this result, we combined the runs into a 
single data set whose rate as a function of LSP is shown in 
Fig.Ill The mean rate is 2.36 ±0.06 events per 10^^ POT 
and the uncertainties shown in the figure are statistical. 

We searched for a sidereal signal by looking for ex- 
cess power in the FFT of the data in Fig. [T] at the 
frequency corresponding to exactly 1 sidereal day. We 
used two statistics in our search, pi — y^Sf + Cf and 
P2 — \/ S2 + C|, where is the power returned by 
the FFT for sm{ujq)TQ), Cf is the power returned for 
cos (ojqTq), and 5*1 and C| are the analogous powers for 
the second harmonics. We used the quadratic sum of 
powers to minimize the effect of the arbitrary choice of 
zero point in phase at O'* LST. Table HIl gives the pi and 
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FIG. 1: Event rate as a function of LSP for the total data set. 



P2 values returned by the FFT for the total data set. 



TABLE II: Results for the pi and p2 statistics from an FFT 
of the data in Fig. [T] The third column gives the probability, 
Vt, the measured value is due to a statistical fluctuation. 



statistic 


p{FFT) 




Pi 

P2 


1.09 
1.13 


0.26 
0.24 



We determined the significance of our measurements 
of pi and p2 by simulating 10"* experiments without a 
sidereal signal. To construct these experiments we used 
the data themselves by randomizing the LSP of each CC 
event 10^ times and placing each instance into a different 
phase histogram. We next randomized the LSP of each 
spill 10"* times and placed the POT for each instance into 
another set of histograms. For the POT histograms we 
drew the phases randomly from the LSP histogram of the 
start times for all spills. Dividing an event histogram by 
a POT histogram produced one simulated experiment. 
The randomization of both the spill and event LSP re- 
moved any potential sidereal variation from the data. 

We performed the FFT on the simulated experiments 
and computed the pi and P2 statistics for each. The re- 
sulting distributions of pi and p2 are nearly identical as 
shown in Fig.H The sofid line at p{FFT) = 2.26 divides 
the P2 histogram at the point where 99.7% of the en- 
tries have lower values of the statistic. Consequently we 
took p{FFT) = 2.26 as the 99.7% confidence limit that a 
measured p{FFT) for either harmonic indicates the dis- 
tribution had no sidereal signal. That is, we adopted 
p{FFT) > 2.26 as our signal detection threshold. 

Based on this threshold, the pi and p2 statistics in Ta- 
ble |TT] show no evidence for a sidereal signal. Thus, the 
normalized neutrino event rate exhibits no statistically 
significant variation that depends on the direction of the 
earth-based neutrino beam in the sun-centered inertial 
frame. In the context of the SME, this result is incon- 
sistent with the detection of LV and CPTV. The third 
column of Table HIl gives the probability, Vf, that the 
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FIG. 2: Distributions of the pi and p2 statistics from the 
FFT analysis of 10'' simulated experiments without a sidereal 
signal. The adopted signal detection limit is p(FFT) = 2.26. 

harmonic powers we found were due to statistical fluctu- 
ations. Vf is the probability of drawing a value of pi or 
P2 from the parent distribution in Fig. [2] at least as large 
as found in the data. 

We next determined the minimum detectable sidereal 
signal strength we could have found in this experiment 
by injecting a sidereal signal of the form Asuy{ojqTq)), 
where A is a fraction of the mean event rate, into a new 
set of 10'' simulated experiments and repeating the FFT 
analysis. We found 68% of the experiments gave pi > 
2.26 for A = 9%. These simulations indicated there was 
a 68% probability of detecting a sidereal modulation if 
the signal amplitude had been 9% of the mean event rate. 

We tested the sensitivity of our results to several 
sources of systematic uncertainty. First we confirmed 
through simulated experiments that the analysis was in- 
sensitive to our choice of zero point in phase, O'' LST, 
due to the definition of pi and p2- Next we tested sources 
that could introduce false sidereal signals into the data 
and mask an LV detection. Degradation of the NuMI 
target caused secular drifts of ~5% in the neutrino pro- 
duction rate on a time scale of ~6 months. Doubling 
this trend introduced no detectable signal for either sec- 
ular decreases or increases. The known ±1.0% uncer- 
tainty in the number of POT per spill Q was too small 
to affect this analysis. Because of the non-uniformity of 
the data taking throughout the solar year, diurnal ef- 
fects, like temperature variations on the POT counting 
devices, could have introduced a false signal. However, 
systematic differences between the day and night event 
rates were smaller than the statistical errors in the rates 
themselves and could not introduce a false signal. Atmo- 
spheric effects could also have imprinted a sidereal signal 
on the data if there were a solar diurnal modulation in 
the event rate that beats with a yearly modulation [l3| . 
Using methods described in [31 , we found that this false 
sidereal signal is < 0.2% of the mean event rate and well 
below the detection threshold. 



Since we found no sidereal signal, we determined upper 
limits on the (0^)^^. and (cl)^:^ coefficients that describe 
LV and CPTV in the SME using the standard MINOS 
Monte Carlo simulation [3]. Neutrinos are simulated by 
modeling the NuMI beam line, including the hadron pro- 
duction and subsequent propagation through the focus- 
ing elements, hadron decay in the 675 m decay pipe, and 
the probability that any neutrinos produced will inter- 
sect the FD 735 km away. These neutrinos along with 
weights determined by decay kinematics, are used in the 
detailed simulation of the FD. 

To find the limits, we chose lAml^l = 2.43 x 10"^ eV^ 
and sin^(26'23) = 1, the values measured by MINOS 
Our tests show that changing these values within the 
allowed uncertainty does not alter the limits we found. 

We determined the limits for each SME coefficient in- 
dividually. We constructed a set of experiments in which 
one coefficient was set to be small but non-zero and the 
remaining coefficients were set to zero. We simulated 
a high statistics event histogram by picking events with 
a random sidereal phase drawn from the distribution of 
start times for the data spills and weighted these sim- 
ulated events by both their survival probability and a 
factor to account for the different exposures between the 
data and the simulation. Simultaneously we simulated a 
spill histogram by entering the average number of POT 
required to produce one event in the FD, as determined 
from the data, at the sidereal phase of each simulated 
event. The division of these two histograms resulted in 
the LSP histogram we used to compute the pi and p2 
statistics. We then increased the magnitude of the non- 
zero SME coefficient and repeated the process until either 
Pi or p2 was greater than the 2.26 detection threshold. 
To reduce fluctuations we computed the limit 100 times 
and averag ed the results. Table Hill gives the mean magni- 
tude of the coefficient required to produce a signal above 
threshold. This procedure could miss fortuitous cancel- 
lations of SME coefficients. We did not consider these 
cases. We cross-checked these limits by simulating 750 
low statistics experiments for each coefficient limit given 
in Table Hm Each experiment had the same total number 
of neutrinos as the data. The distributions of the pi and 
P2 statistics for all the experiments were used to deter- 
mine the confidence level at which the measured values 
in Table HI] are excluded by each limit. The exclusion 
using this method is > 99.7% C.L. for all coefficients. 

In summary we found no evidence for sidereal varia- 
tions in the neutrino rate in the MINOS FD. This result, 
when framed in the SME 0, @] , leads to the conclusion 
that we have detected no evidence for the violation of 
Lorentz or CPT invariance described by this framework 
for neutrinos traveling over the 735 km baseline from 
their production in the NuMI beam to the MINOS FD. 
The limits on the SME coefficients in Table Hill for the 
FD that come from this null result improve the limits we 
found for the ND by factors of order 20 - 510 6]. This 
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TABLE III: 99.7% C.L. limits on SME coefficients for 

Vt', {O'Lj'^T have units [GeV]; {cl)%t are unitless. The columns 

labeled X show the improvement from the near detector limits. 



Cocff. 


Limit 




I 


Cocff. 


Limit 




I 




5.9 X 10" 


-23 


510 




6.1 X 10" 


23 


490 




0.5 X 10" 


23 


20 




0.5 X 10" 


23 


20 




2.5 X 10" 


23 


220 




2.4 X 10" 


23 


230 




1.2 X 10" 


23 


230 




0.7 X 10" 


23 


170 




0.7 X 10" 


23 


190 











impi'ovcimcnt is due to the different behavior of the oscil- 
lation probability in the short and long baseline approxi- 
mations coupled with the significantly increased baseline 
to the FD. These improvements more than offset the sig- 
nificant decrease in statistics in the FD. They are the first 
limits to be determined for the neutrino sector in which 
LV and CPTV are assumed to be a perturbation on the 
conventional neutrino mass oscillations. 
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